Quantum coherence and carriers mobility in organic semiconductors 
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We present a model of charge transport in organic molecular semiconductors based on the effects of 
lattice fluctuations on the quantum coherence of the electronic state of the charge carrier. Thermal 
intermolecular phonons and librations tend to localize pure coherent states and to assist the motion 
of less coherent ones. Decoherence is thus the primary mechanism by which conduction occurs. It 
is driven by the coupling of the carrier to the molecular lattice through polarization and transfer 
integral fluctuations as described by the hamiltonian of Gosar and Choiii Localization effects in the 
quantum coherent regime are modeled via the Anderson hamiltonian with correlated diagonal and 
non-diagonal disordei"^ leading to the determination of the carrier localization length. This length 
defines the coherent extension of the ground state and determines, in turn, the diffusion range in 
the incoherent regime and thus the mobility. The transfer integral disorder of Troisi and Orlandi^ 
can also be incorporated. This model, based on the idea of decoherence, allowed us to predict the 
value and temperature dependence of the carrier mobility in prototypical organic semiconductors 
that are in qualitative accord with experiments. 
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I. INTRODUCTION 

Molecular electronics is a field that is rapidly gaining 
importance because of its potential in producing a new 
breed of plastic organic devices. A physical quantity that 
is critical to their operation is the charge carrier mobil- 
ity in the organic molecular semiconductor that forms the 
active layer in these devices. It is therefore not surprising 
that along with this resurgence of interest in molecular 
electronics in recent years, charge transport mechanisms 
in organic semiconductors are once again at the forefront 
condensed-matter physics research. Consequently, im- 
portant debates in this field have been rekindled recently 
to which we would like to contribute by the present work. 

Experimental observations of " band-like" charge trans- 
port in single crystalline organic molecular semiconduc- 
tor have previously been observed in stilbene^, perylene*^, 
acenes^ or rubrenCf^iS, Hole mobilities of 1 to 10 cm^/V.s 
at room temperature were reported. In many cases, 
these mobilities were shown to decrease with increas- 
ing temperatures. At low temperatures, values of a 
few hundred cm^/V.s were even seen in a "time-of- 
flight" measurement 1^ Earlier theoretical explanations of 
these results invoked rigid band models to account for 
the power-law temperature dependence of the mobility.— 
Later on, the electron-phonon interaction became widely 
recognized as the key factor that determines the mobility. 



leading to the development of several polaronic models 
of charge transport in these materials ] ^"'"'^"'^i-'^^'-'^'^i-'^^'-'^^i"'^^ 

The first proposed polaronic models have emphasized 
the influence of the molecular character of organic molec- 
ular crystals on charge transport. One approach calcu- 
lated the reorganization energies on a single molecule and 
used the Marcus theory of charge transfer-iSiiiii^iH to de- 
termine the mobility. Another— li^ii^ stressed the impor- 
tance of low energy phonons and librations and derived 
an average mobility from the Kubo formula. In their re- 
cent works, Troisi and Orlandi criticized both these ap- 
proaches based on polaron theories. They have suggested 
that the formation of a small polaron with the charge 
localized on a single molecule is unlikel}*^^ because the 
nuclear reorganization energy in these crystals is compa- 
rable to the average intermolecular charge transfer while 
the more general approach based on intermolecular vibra- 
tions suffers some problems of averaging;^ Subsequently, 
these models were extended to describe the effects of ther- 
mal disorder and thermal fluctuations on charge trans- 
port. In this case, low energy phonons and librations are 
able both to localize the charge and to drive its diffusion 
in the lattice. In a previous work)^ we implemented the 
idea of phonons acting as a source of disorder and lo- 
calizing the charge using an Anderson hamiltonian with 
the carrier coupled to low energy intermolecular phonons 
and librations through the fluctuations of the polariza- 
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tion energy. On the other hand, the model of Troisi and 
Orlandi focused on the effects of fluctuations of the trans- 
fer integral which are considered to be large enough to 
localize and then move the carrier in the latticed We note 
that these theoretical considerations^i^ have been intro- 
duced in the pioneering work of Gosar and Choi! on the 
determination of the mobility of an excess charge in a 
molecular crystal. 

The present work focuses on quantum coherence ef- 
fects on charge transport in organic semiconductors that 
has been ignored so far because of the technical diffi- 
culty of treating decoherencc in a system with strong 
electron-phonon interactions. We calculate the carrier 
localization length resulting from the fluctuations of the 
electronic polarization energy. The fluctuations of the 
transfer integrals treated in Ref. can also be included 
without problem. This localization length allowed us to 
determine the carrier diffusion coefficient leading to the 
prediction of the mobility of charge carrier in organic 
molecular crystals. We find that these results can explain 
the main trends of time of flight mobility measurements 
performed in pure organic semiconductor single crystals. 



II. THE GENERAL FRAMEWORK OF THE 
CONDUCTION MODEL 

Low energy phonons and librations play an ambiva- 
lent role in limiting charge carrier transport in molecular 
semiconductors. Depending on the quantum coherence 
of the electronic state of the carrier, thermal phonons 
will tend to localize pure coherent states or assist the 
motion of less coherent ones. As long as the carrier elec- 
tronic states keep their quantum coherence, they are es- 
sentially localized by thermal disorder within a localiza- 
tion length i, which defines the coherent extension of 
the ground state. At longer time scales, thermal fluctu- 
ations cause decoherence and assist the diffusive motion 
of the carrier. This process is driven by the coupling 
of the carrier to the molecular lattice, through polar- 
ization fluctuations and transfer integral fluctuations as 
described by Gosar and Choii The localization length L 
determines the diffusion range and thus the mobility. To 
determine this quantity in the quantum coherent regime, 
we first introduce a transfer matrix formalism in two di- 
mensions depicting the quantum interference processes 
in the plane of high conduction of these organic semi- 
conductors. The hamiltonian for this quantum model is 
the Anderson hamiltoniani^ with correlated diagonal and 
non-diagonal disorder established in Ref. @, and which is 
extended in the present work by adding the purely diag- 
onal transfer integral disorder computed in Ref. |3|. The 
role of spatial and energetic correlations will be stud- 
ied with particular emphasis for the case of polarization 
fluctuations. For the sake of concreteness, we shall ap- 
ply the results derived here to widely studied molecular 
single crystals such as the acenes or rubrene. 

The lattice dynamics in these materials have been in- 



vestigated extensively. Inelastic neutron scattering in 
naphtalene by Natkaniec et alii^ revealed twelve inter- 
molecular phonon branches in crystalline acenes. Struc- 
tural refinements using a molecular dynamics model 
yielded a satisfactory fit to the measured phonon disper- 
sion curves. More recently, some of these modes along 
with their relative electron-phonon coupling constants 
were obtained by a density functional theory scheme 
(DFT/LDA) applied to the acene series from naphtal- 
ene to tetraceneji^ Both these calculations converge on 
the existence of a vibrational band centered near 50 cm~^ 
with 3 acoustic modes, 3 optical modes, and 6 librations. 
The directional average of the root-mean-square ampli- 
tudes of the translational vibrations of the molecules is 
0.17A at room temperature in anthracene as obtained 
from the X rays measurements of Gruickshanki^ This 
work also showed a typical librational amplitude of 3 or 
4 degrees at room temperature. The largest hole transfer 
integral between adjacent sites in pentacene, the proto- 
type acene compound, is about 1000 cm~| (0.12 eV) as 
determined independently in Refs. ITtI and [2ll. Thus, the 
thermal motion can essentially be considered as static 
within the hole residence time of the order of 0.1 ps. 
More details on time scales are given in Appendix [X] 
This is enough time to build a large quantum coherence 
into the system based on quantum interference in the An- 
derson weak localization regime. For times longer than 
1 ps, the coupling of the extra charge with the inter- 
molecular phonons reservoir now becomes the source for 
decohcrenccj^ In this case, we treat the charge motion 
within an adiabatic classical diffusion approximation, as 
the charge follows the motion of the low energy phonon 
wavepacket. The characteristic length of this diffusion 
is just the localization length determined in the short 
time scale coherent regime. Intramolecular vibrations 
play a negligible role in the decoherence. In fact, ac- 
cording to Ref. [ll|- there is essentially one mode that 
contribute to the reorganization energy of the positively 
ionized molecule. In pentacene, this corresponds to vi- 
brational modes at 1340 cm~^. This mode is too fast to 
localize the charge. Its main effect is to renormalize the 
transfer integral as outlined in our calculation presented 
previously in Appendix C of Ref. [H. That calculation 
yielded a further reduction of the bare transfer integral 
by a factor of 0.75. 



III. THE CALCULATIONS OF THE 
LOCALIZATION LENGTHS 

A. The hamiltonian 

On each site of a real molecular crystal at finite tem- 
perature, the molecular positions and angles fluctuate 
with respect to their values in the perfect crystal. The 
carrier motion in the lattice is determined by its inter- 
action with these fluctuations. In order to construct 
the hamiltonian that properly accounts for these interac- 
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tions, a clear distinction must be made between fast and 
slow interactions with respect to the transfer time h/J^, 
where J° is the larger transfer integral in the conducting 
plane of the molecular semiconductor. This classifica- 
tion of the interactions according to their time scale is 
described in Appendix |^ we show that the interaction 
with the electronic polarization and with intramolecular 
vibrations are fast, while the interaction with intcrmolec- 
ular phonons and librations is slow. Fast interactions can 
be averaged; they just renormalize the parameters of the 
hamiltoniani^ On the contrary, slow or static interac- 
tions have to be included explicitly in the hamiltonian. 

This is the case for intermolecular phonons, which as 
we mentioned earlier, form a band around 50cm^^ in 
acenes. These displacement fields on each crystal site 
can be considered as uncorrclatcd because of the nearly 
degenerate twelve phonon modes of acenes : many low- 
frequency phonons can also be treated as independent 
random displacements Srj and S9j on each site j . These 
lattice fluctuations as shown below, induce, in turn, site 
energy fluctuations Sj and transfer integral fluctuations 
S'j jj^h ■ These random variables are characterized by 
mean square values 5 and 5' respectively. We can now 
start the construction of the hamiltonian. 

In a previous paper—, the long range Coulomb hamil- 
tonian describing the Coulomb polarization induced by a 
charge in the molecular lattice is mapped onto a short 
range tight binding hamiltonian containing the polar- 
ization energy Ep(n) on each molecular site n and the 
renormalizcd transfer integral Jn,n+h which couples two 
adjacent sites n and n + h. 



cloud and the intramolecular phonon cloud which dress 
the charge. The reduction factors of the bare integral J 
have been calculated in Ref. [H for pentacene to be 0.79 
for polarization and 0.75 for intramolecular phonons. 

The situation changes greatly when intermolecular 
thermal disorder enters the system. Thermal disorder 
on each site is parameterized by six Gaussian random 
variables attributed to each translational and librational 
degree of freedom of the molecules. 

We designate the position fluctuations 6Rj ~ fj — ff' 
and the angular fluctuations arc represented by 69 j — 

9j — 9j . Following Ref. [3, the polarization energy in the 
presence of thermal disorder can be written at each site 
n as 



Ep{n) ^ E^pin) + Sr^ (2) 
where (5„ = Y.j^ni^^pi^)/^'^!) i^^li ^ ^'^n) + 
{dEp{n)/d9j) {69 j — 59n)- The renormalization factor of 
the transfer integral at each site due to polarization ef- 
fects is given by^ 



J n.n 



+h 



n.n+h 



exp ■ 



5n + 6n+h 



(3) 



where A, an energy of the order 1 eV, defines the im- 
portance of renormalizationi^ It is also possible to in- 
troduce in the hamiltonian the fluctuations of the bare 
transfer integrals, Jn.n+h, following the model of Troisi 
and Orlandii as 



H 



^ Ep{n)\n) {n\ - ^ Jn^n+h\n){n + h\ (1) 



n.n-\-h 



Jn,n+h{rn+h — r9n+h ^ 9n) 



jO 

^n,n-\-h 



n,n-\-h 



(4) 



In the case of a perfect crystal, the polarization energy 
is uniform and shifts the ground state uniformly by about 
1 eV as observed experimentally in acenes i^*^'^^ Conse- 
quently, the bare bandwith is significantly narrowed, re- 
gardless of the temperature, both by the polarization 



where (5^ 



d{6rn+h+Srn) + iiS9n+h-S9n). Then 
hamiltonian [T] becomes an Anderson hamiltonian with 
correlated diagonal and non diagonal disorder and in- 
cludes both polarization fluctuations, and transfer 
integral fluctuations, S'^ 



H 



J2iEUn) + 5n)\n){n\ 



E 

n,n-\-h 



^n,n+h)cxp 



5n + <5,i+h 



n){n + h\ 



(5) 



r 



In principle, the integral energy fluctuations 5'^ 
and the polarization fluctuations (5„ are correlated be- 
cause they are derived from the same thermal displace- 
ment field ((5r„,i50„). It is possible that the complexity 
of the displacement fields dilutes the effects of these cor- 
relations. Troisi and Orlandi have studied the effects of 
5'^ ri+ft, alone. Here we shall focus on the polarization fluc- 
tuations 5m which appear both in the diagonal and the 



nondiagonal part of the Anderson hamiltonian. Thus in 
the following we shall set 5'^ = in the hamiltonian 
of Eq. [H We are left with the random variable 6n the 
distribution of which has a root mean square value 5. 

Polarization fluctuations effects are long ranged. In- 
deed the polarization cloud extends over many molecules 
as the contribution of the induced dipoles decrease like 
the reciprocal distance to the carrier (l/r). It is thus im- 



4 



TABLE I: Thermal energetic disorder S resulting from the 
translational, librational and both translational and libra- 
tional lattice fluctuations is the root mean square value of the 
disorder distribution. The renormalization factor A indepen- 
dent of the type of disorder (see Eq. [5]) is related essentially 
to the HOMO-LUMO gap^ 



Disorder 


S (meV) 


A (eV) 


0.1 A 


42.1 


~ 0.4 


3 degrees 


9.4 


~ 0.4 


3 degrees + 0.1 A 


44.9 


~ 0.4 



portant in our calculations to explore the amplitude of 
the correlations. This has been done in Appendix [Bl The 
result is that spatial correlations due to these disordered 
dipoles are much shorter range. Thus the correlations be- 
tween the polarization energy of two adjacent sites can be 
neglected. Nevertheless, diagonal and nondiagonal disor- 
der are correlated through the renormalization factor A. 
The disorder parameters which constitute the starting 
point of the transfer matrix calculation are presented in 
tablelH The root mean square values 5 of the polarization 
fluctuations have been calculated as follows. For a given 
distribution of lattice displacements that we know from 
the experimental work of Cruickshank2£^ the polariza- 
tion energy has been calculated according to our previous 
papers.i^^ Then by varying the number of samples, the 
distribution of (5„ can be reconstructed. We have checked 
that it is gaussian with a root mean square energy S. 



in favor of the backwards one. Such a problem cannot 
be treated by averaging the energy distribution because 
the coherence of the quantum state, expressed by the 
phases of the wave function at each molecular site, should 
be preserved. The only transport model which is com- 
pletely quantum is the transfer matrix formalis m^^i^^ , 
which determines, in amplitude and phase, the carrier 
transmission and reflexion coefficients in a disordered lat- 
tice. These coefficients cannot be directly averaged. Only 
the Lyapunov of the distribution obeys a central limit 
theorem and can be averaged i^^i^^'^^ We used the trans- 
fer matrix formalism to calculate the localization length 
of the carrier. In two dimensions this method can be 
applied to a long strip. 

We consider a 2D lattice array system of size N x M 
sites as illustrated in Figure [TJ Each site (n, m) corre- 
sponds to one molecule and is characterized by an en- 

N 



M 



FIG. 1: Two-dimensional strip of width A4 and length A'^. 



ergy (5„,„ from a gaussian distribution. The ground state 
wavef unction lib) 



B. The localization length 

The hamiltonian of the carrier coupled to slow lattice 
fluctuations of Eq. [5] is of the Anderson type (with cor- 
related diagonal and nondiagonal disorder). It is well 
known that such types of hamiltonian lead to Ander- 
son localization: the carrier quantum state results from 
the interference process between the wavelets scattered 
at each site; this process weakens the forwards scattering 



IV") = ^a„^m\n,m) (6) 

is solution of the Schrodinger equation 

H|V^) = E\yj) (7) 

We impose periodic boundary conditions and only take 
into account next neighbor interactions. Then equation 
[7] can be written explicitly as. 



{n,m\{H - E)\i') = 



(8) 



J" 



r,m^ 

n.rn "n— 1,' 



Using An = (an,i, . • . , an,M) then the former equation 
becomes, 

An+l = MnAn + Mn-lAn-1 (9) 

with M„ and M„ given in Appendix [Cl This equation 
can be recast into the transfer matrix formalism as. 



f A„+i \ _ f Mn+l Mn \ { An \ _rr. ( An \ 
\ An I n )\ An-l J ~ n ^"-1 / 

(10) 

We define —Y\^ T„ then 
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An+l 
An 



(11) 



For the limiting case of n going to infinity, t„ exposes 
the asynitotic behaviour of the ground state wavefunc- 
tion According to Oseledec's theoren*^^^ the foUow- 
ing hmit yields the quantity 7^, called i-th Lyapunov ex- 
ponent where Vi is i-th eigenvalue of r„. 



hm = -7i 



n^oo TL 

is well defined so that when n increases , Vi is simply 



(12) 



-n7i 



(13) 



Now, 7i > describes the exponential rate of decrease 
of the components of the ground state |^). Then, any 
solution of the Schrodinger equation will decrease faster 
than the minimum rate given by 



7 = min(7,;) 



(14) 



L = 1/7 is called localization length and represents a 
measure of the extension of the ground state wavefunc- 
tion. 

To calculate the localizaton length, we sorted 2N ran- 
dom vectors and calculated their evolution via the trans- 
fer matrix process. Actually, if we start from a vector Vq, 
after n steps, its image Vn is dominated by the largest 
Lyapunov exponent. To obtain the smallest one, that 
corresponds to the localization length, we used an or- 
thogonal normalization^ process applied to 2N random 
vectors in order to get rid of the contribution of the 27V— 1 
first Lyapunov exponents that screen the one we look af- 
ter. 

These Lyapunov exponents were calculated for arrays 
of finite width M . We have to use a finite-size rescaling 
to extract the actual localization length L from the lo- 
calization length L{M) for finite values of M. According 
to Ref. Hll, as the diagonal parameter 3/ J is lower than 
4, we identify in our 2D system the localization length L 
of the infinitely wide array by taking the limit 



M 



lim L{M) 



(15) 



where L{M) is the localization length obtained for an ar- 
ray of width M. Figure [2] shows the values of L{M)/M 
versus 1/M for different energy positions in the band- 
width calculated for A = 0.4 as established in Ref. for 
an angular disorder. As 1/M goes to 0, we can fit, as 
represented on figure [21 L{M)/M (that also goes to 0) 
with a second order polynomial a/M + b/AP to find the 
expected localization length L equal to the parameter a. 



A complete study of such a 2D system was recently 
carried out by Unge et Stafstrom. The results of Ref. [2^ 
were partly used to validate our procedure of calculation 
of the Lyapunov exponents. 
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FIG. 2: Localization length in a strip of finite width L{M). 
Each curve corresponds to a diflferent energy in the band and 
a disorder of 0.5 J. When 1/M tends to 0, the locahzation 
length for a given energy corresponds to the asymptotic slope 
of the corresponding curve. 

Since we are interested in charge transport, we explore 
the behavior of the localization length near the band 
edge. Figure [3] represents the results for an energy range 
between — 3.95J and — 4.05J with the band-edge of the 
perfect 2D lattice being equal to — 4J. 



IV. THE EVALUATION OF THE MOBILITY 

At time scales of a fraction of picosecond or less, inter- 
molecular phonons can be viewed as static disorder; then, 
at this time scale, quantum coherence induces carrier lo- 
calization. A carrier placed in such a system extends 
coherently over a characteristic length L. The question 
is how can this electron that is localized by thermal dis- 
order move to a new position? If all types of intramolec- 
ular phonons wether acoustic, optical and librational, are 
able to localize the charge efficiently, they cannot drive 
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FIG. 3; Localization lengtli close to the hole band-edge for 
different values of the diagonal disorder. When the disorder 
increases, the localization length decreases. 



the motion of the charge with the same efficiency. 

Dispersionless optical phonons are not efficient in mov- 
ing charges. Emin was the hrst to observe this facti^ 
which was established more rigorously in Ref. [s^ Fur- 
thermore wc consider the charge motion to be adiabatic, 
because transfer integrals are at least one order of mag- 
nitude larger than intermolecular phonon energies. The 
most efficient process for moving the extra-charge is thus 
to use an acoustic phonon moving in the conduction 
plane. Then, during the lifetime of the localized excita- 
tion, the charge follows the phonon packet adiabatically. 
It moves along with the acoustic phonon with a speed 
that is basically equal to the sound velocity Vs- Values of 
the sound velocity in naphtalene are given in Ref. H and 
can also be deduced from the phonon dispersion curves 
of Ref. [l^. We have taken a value of Vg = 3.3km.s~^ , 
which we consider a good order of magnitude for acenes. 

Thus the decoherence from the Anderson localized 
state to the classical diffusion in the plane, leads to a 
mobility 
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(16) 



The Einstein's relation used here is valid for a non- 
degenerated hole gaz. This condition is achieved in time 
of flight measurements and in field effect transistors. 
From the transfer matrix calculation applied near the 
hole band edge, we find a localization length L/a vary- 
ing as a power law of the polarization energy fluctuation 
(cf Fig. Il 
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(17) 
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FIG. 4: Neper logarithm of the localization length calculated 
at -E = — 4J versus the Neper logarithm of the diagonal dis- 
order. The linear relation leads to a power-law behaviour as 
in experiments. The slope is about —1.4. 



where a = 5 nm is the lattice spacing parameter and 
J ~ 58 mcV the renormalized transfer integral. 

If we combine the effect of a Gaussian translational 
disorder with a root mean square amplitude of 0.1 A and 
a gaussian librational disorder of 3 degrees around each 
axis, we get from table |I] a resulting energetic disorder of 
45meV. From Figure[31 we deduce a localization length of 
11a. Using Eq. we get a mobility of 1.8 cm^.V^^.s^^ 

Furthermore, because the polarization coupling is lin- 
ear in the different molecular degrees of freedom (Eq. 
[2]), the fluctuation S is proportional to (/cbT)^/^ at high 
temperatures. Then, we find that the mobility varies as 
a power law of the temperature with exponent a 2± —1.7. 

It is important to have an idea of the temperature 
range of validity of the power laws that results from 
our theory. They are based on a diffusion-like view 
of the transport. The disorder characterized by (5 is a 
site disorder, whereas the localized wavcfunction spans 
multiple sites and is therefore characterized by a disor- 
der that is lower by a factor of (a/L)^. The condition 
(5(a/L)^ < ksT guarantees the diffusive aspect of the 
motion; otherwise, for cases in which the effective disor- 
der exceeds fcsT, a hopping process appears. With the 
values of Table U concerning thermal disorder and polar- 
ization fluctuations this condition yields T > 50 K. 



V. COMPARISON WITH EXPERIMENT 

The present theoretical work is applicable to the eval- 
uation of the intrinsic mobility of a carrier, which prop- 
agates through a single crystalline molecular semicon- 
ductor such as pentacene or rubrene. In particular, it 
is not directly applicable to the channel of a molecular 
transistor with an oxide gate. In these transistors, sur- 
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face Frohlich polarons have been shown to play a major 
role.— Moreover, the field-effect mobility is affected by 
the presence of traps, particularly at the gate interface. 
Charge carrier traps dominate transport in thin film tran- 
sistors and often in single crystalline transistors These 
traps have not been considered here. 

The temperature range of interest is between 50 and 
400 K where the mobility decreases with increasing tem- 
perature with a power law, T^". The relevant exper- 
iments were performed on large ultrapurified molecular 
crystals by using the time-of-flight technique. The expo- 
nent a was found to be 2.9 for holes and 1.4 for electrons 
in naphtalene,^^ 1.5 for holes and 1.26 for electrons in 
anthracenes^ and 1.87 for electrons in perylenci^ More 
recently the measured values of a in sublimation grown 
perylene single crystals were found to be 2.8 for electrons 
and close to zero for holes.— In biphenyl single crystals, 
the electron mobility exponent is 1.18 for both electrons 
and holesi^ In phenenthrcne, the hole mobility varies with 
the exponent 1.8 and the electron mobility with the ex- 
ponent 0.95.— In tetracene single crystals, an exponent a 
of 2 was reported for holes^, while for air-gap transistor 
in rubrene a value of 1.45 can be deduced above 240 K 
from the curve of Fig. 3 in Ref. [ssl . 

All the experiments cited above demonstrate that in ul- 
trapure single crystals of molecular semiconductors, the 
mobility decreases as a power law with increasing tem- 
perature. The precise value of the exponent depends on 
the details of the crystal structure, on the transfer in- 
tegrals, on the polarizabilities of the molecules, and on 
the interaction of the extra charge with the lattice. The 
exponent a varies essentially between 1 and 2 depending 
on these factors. 

In this sense our theoretical model, which contains no 
adjustable parameter and describes a prototypical crys- 
talline molecular semiconductor is in excellent agreement 
with experiment. The model can be easily adapted to 
calculate a in many systems by entering the relevant 
material-specific parameters. The geometry of the lat- 
tice can also be chosen according to more precise crystal- 
lographic data (in the present analysis, the conducting 
plane of acenes is assumed to be a square molecular lat- 
tice). The values of the relevant transfer integrals and 
their number can also be varied (here we have consid- 
ered a single transfer integral). In general, these val- 
ues are different for electrons and holes. The calculation 
of thermal disorder at a given temperature can also be 
varied according to the polarizabilities of the individual 
molecules and the characteristics of the coupling to the 
lattice. Moreover, the presence of defects and impurities 
acting as traps in the actual crystals can always change 
the temperature dependence of the mobility, even for 
temperatures very outside of the hopping regime. Only 
intrinsic effects have been included in the present calcu- 
lation. 

It is worth summarizing here the different parameters 
that we have used in the model and recalling their origin. 
Most of these values concern pentacene, but in the cases 



where the values does not exist for this compound, we can 
also infer the value from experiments on other acenes. 

- bare transfer integral : ~ 100 meV from Refs. [2ll 
and 113 

- effective transfer integral renormalized twice accord- 
ing to Ref. [2I, J = 58 meV 

- relevant intramolecular phonon frequency j^i 1340 
cm~^ 

- typical intermolccular phonon frequencies f^^*^ 50 
cm~^ 

- sound velocity in the plane jiii^ 3,3 km/s 

- amplitude of the translational vibrations at room 
temperature)^ 0,17A 

- amplitude of the librations at room temperature)^ 
'-^ 3 degrees 

- polarization energy of a carrier Ep ~ —1.5 eV, cal- 
culated in Ref. [H. 

Special cases have also been observed where the mobil- 
ity is relatively constant over a large temperature range. 
In field effect transistors this is now attributed to the 
existence of Frohlich polarons i^^i, In the present form our 
model cannot account for such a low value of a (close to 
zero). 

VI. COMPARISON WITH RECENT THEORIES 

The present theory depicts the behavior of a carrier 
in a random medium resulting from lattice fluctuations. 
It is important to compare it to other existing models 
that pursue the same purpose by using other scenarios: 
polarons, reorganization energies, etc. 

In our model, lattice fluctuations are coupled to the 
charge carrier through fluctuations of the Coulomb po- 
larization energy— The idea is to consider possible dy- 
namic localization processes induced by these Coulomb 
fluctuations acting on the carrier as a random fleld. This 
means that at short time scales, the motion of the carrier 
is coherent and can be treated quantum-mechanically us- 
ing the transfer matrix formalism. At this stage we avoid 
any averaging of the coupling energies. We keep quan- 
tum interference and Anderson localization, which arc 
relevant to this problem. A study of Anderson localiza- 
tion in molecular semiconductors using the same trans- 
fer matrix formalism was recently published by Unge and 
Stafstromi^ The hamiltonian they use is the same as Eq. 
[5] and comes from Ref. S Our results are consistent with 
theirs, but the authors do not follow the same analysis 
or present a transport model. 

We consider the localization process to be dynamical. 
At time scales longer than 1 ps (see Appendix]^ decoher- 
ence occurs. Then the localized carrier "surfs" adiabat- 
ically on the acoustic phonon waves: in fact, the carrier 
localizes in a certain landscape but the landscape moves 
slowly and the carrier is forced to follow this motion. 

In some respect, this picture resembles the scenario 
proposed by Troisi and Orlandi^ to achieve the same pur- 
pose. Our model is 2D instead of theirs which is ID 
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and we use quantum transport instead of a semiclassical 
simulation. Both models insist on the fact that lattice 
fluctuations can both localize the carrier and drive its 
diffusion in the lattice. 

It is also important at this stage to compare our theory 
with well established recent models based on reorganiza- 
tion energies and polarons. 

In these last years, Marcus theory has been extensively 
used to deduce charge transport parameters in molecu- 
lar organic semiconductors. The model has been estab- 
lished half a century ago to predict charge transfer rates 
between a reactant and a product in a donor-acceptor 
reaction. More recently it has been widely used to de- 
termine small polaron hopping rates in oxides such as 
chromia, iron oxide, anatase and rutile,— and in molec- 
ular organic scmiconductorsJ^iiiii^ii^ The advantages of 
this model are that the transport parameters are consid- 
ered to depend only on a pair of adjacent sites and that 
the carrier is coupled to the lattice through optical or 
intramolecular vibration modes only. The result is sim- 
ple and the charge transfer rate kij can be written in the 
semiclassical and nonadiabatic approximation : 

k --t -^ /^^exp ^i^^^l^M!^ (18) 

where Xij is the reorganization energy, AEij = Si — Sj, 
Ei and Ej are the energies of the initial and final states, 
and tij is the transfer integral between the two sites. 

The value of Ay corresponds to the dimer energy dif- 
ference between the situation where the pair is charged 
but the geometric configuration corresponds to a neutral 
pair, and the situation where the pair is charged and is 
in the true geometry. 

It is worth noting that, in this form, Marcus theory is 
strictly equivalent to the Emin-Holstein's earlier model of 
small polaron hopping.— Emin and Holstein introduced 
the concept of polaron binding energy Eb (equal to one 
quarter of the reorganization energy) and the idea of co- 
incidence, which is equivalent to the idea of a transition 
state in Marcus theory. 

Due to the fact that Marcus theory always yields a 
transfer rate that increases with temperature, it cannot 
adequately describe the charge mobility in the bulk of 
single crystals where the temperature behavior is just 
the opposite. Even in the hopping regime, relation [TH] 
has many restrictions. For example, this equation is only 
applicable when A/4 is much larger than i^j, a condition 
that is not satisfied for acenes and related compounds.— 
It is also only valid when the relevant intramolecular 
phonon frequency responsible for charge transfer (1340 
cm~^) is much lower than the transfer integral tij. This 
is also not the case in acenes and related compounds.— A 
third restriction is that small polaron hopping in a crys- 
tal must obey selection rules.— The optical modes in- 
volved in the small polaron formation cannot lead to car- 
rier motion in the lattice. Acoustic phonons are manda- 
tory for this hopping process to occur. This solid state 



effect is not accounted for in Marcus theory. 

Finally relation [18] includes only short range interac- 
tions which are usually deduced from ab-initio software 
packages. As emphasized in reference "the polar- 
ization effect in these systems is largely electrostatic in 
nature and can change dramatically upon transition from 
a dimer to an extended system." 

In order to overcome the problem of nonlocality and 
to introduce the low energy acoustic phonons that are 
of paramount importance in charge transfer processes in 
molecular semiconductors, Hannewald and Bobbert have 
worked on the basis of a Peierls-Holstein hamiltonian. 
Both the carrier and the phonons are treated quantum 
mechanically. The role of low energy phonons is to mod- 
ulate the carrier energy on each site and the transfer 
integrals between adjacent sites. At finite temperature 
this modulation due to a large number of phonons modes 
could induce Anderson-like localization on the carrier, 
especially in anisotopical electron systems like acenes. 
However, as already observed by Troisi and Orlandi, the 
treatment of Hannewald and Bobbert excludes this pos- 
sibility by averaging the electronic energies at all time 
scales 4 



VII. POLARIZATION ENERGY 
FLUCTUATIONS OR TRANSFER INTEGRAL 
FLUCTUATIONS ? 

In general, thermal energetic disorder in organic molec- 
ular semiconductors comes either from a distribution of 
polarization energies or of transfer integrals or a combina- 
tion of both of these fluctuations. Troisi and Orlandi^iii 
have shown that by an appropriate choice for the val- 
ues of the parameters in their model, transfer integral 
fluctuations alone are able to account for the order of 
magnitude of the mobility and its observed temperature 
dependence. Here, we have also shown that polarization 
fluctuations alone can achieve the same result. 

In fact, the intermolecular potentials contain both 
short range and long range contributions due to the pres- 
ence of an extra charge. It is typical in quantum chem- 
istry and molecular dynamics calculation a'^'^^ that only 
the short range part of these potentials are considered. 
However, long range polarization effects are not negligi- 
ble and must be included in any realistic calculation. 

It is important to emphasize the result of two recent 
experiments which show that polarization contributions 
cannot be avoided and do not represent special cases. 
The measurements of Morpurgo et ali^ have shown un- 
ambiguously that the hole field effect mobilities in sin- 
gle crystals of rubrene or tetracene depend strongly on 
the dielectric permittivity of the gate. We have recently 
shown2^ that such dependence can be understood only 
if polarization effects are taken into account. Our re- 
cent thermopowcr measurements^ on high quality pen- 
tacene films evaporated on different substrates revealed 
a large intrinsic temperature independent, contribution 
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to the Seebeck coefficient of 265/iV/K. This imphcs that 
each carrier transports an intrinsic vibrational entropy of 
3 Boltzmann constants (Sfes)- We have shown quantita- 
tively that this large entropy is associated with local vari- 
ations of the low energy vibration frequencies around a 
carrier by 30% with respect to the bulk. This fact, which 
is not considered in Ref. [3, is easily understood when one 
considers that polarization effects induced by the carrier 
locally change the character of the interaction between 
molecules from Van dc Waals bonds to point-dipole or 
dipole-dipole interactions. They validate the concept of 
polarization fluctuations which was introduced initially 
by Gosar and Choii 

In fact, the role of long range polarization in the motion 
of carriers in acenes and related compounds was adrcssed 
for the flrst time in 1954^1^ and considered seriously in 
the seventies4^ Nowadays it tends once more to be widely 
recognized in the scientific communityj ^^i^°'^^ 



VIII. CONCLUSION 

Following the work of Gosar and Choi^ we have pro- 
posed a model that accounts for the main features of 
transport in semiconducting acenes and related com- 
pounds. Like Troisi and Orlandi)^ we think that charge 
transport is due to lattice fluctuations. An important 
contribution of the present model is the inclusion of quan- 
tum interference effects, which are generally ignored in 
other transport models. 

The present theory exclusively concerns conduction in 
bulk ultra pure crystals. It is also likely to apply to 
conduction in ultrapure crystal transistors made with 
polymer or air-gap gates. Although Frohlich polarons 
were not considered in the present work, we have pro- 
posed that they dominate the transporli^ in cases where 
a molecular single crystal is interfaced with an oxide gate. 
This conclusion is consistent with recent experiments)^ 

Even in single crystals, the presence of electroactive 
defects acting as traps has been demonstrated They 
change the temperature dependence of the mobility. In 
the case where these traps dominate the transport, typ- 
ically in thin-fllm acene transistors, a hopping regime 
sets in and the mobility becomes essentially activated 
with temperature. Even in this regime, the electronic 
polarization induced around trapped and free charges re- 
mains of paramount importance. A work is in progress 
to demonstrate this point theoretically. 



APPENDIX A: TIME SCALES 

The non-interacting band properties of a perfect pen- 
tacene crystal along all crystallographic directions have 
been calculated by Troisi et Orlandi^^ and Cheng et alj^ 
For the bare transfer integral, J°, between molecules 
along the direction of easy propagation in the (a, 6)-planc, 
they obtain — 100 mcV from which one gets hj 



4 X 10""'^ s as the characteristic time for in-plane Bloch- 
wave formation. The corresponding transfer time in 
the perpendicular c-axis direction, /i/J^^, is thirty times 
longer than in the plane. Thus, in the presence of scatter- 
ing, which substantially reduces the Bloch-wavc lifetime, 
the carrier motion is essentially two-dimensional. 

The above considerations allow for the classification 
of the various charge carriers interactions in organic 
semiconductors. For fast interactions with characteristic 
times shorter than /i/J", the charge can be assumed to 
be located on a single molecular site. In pentacene, this 
is the situation encountered during the interaction of the 
carrier with the electronic polarizability of the medium 
or in intramolecular charge-transfer as well as the cou- 
pling with intramolecular carbon stretching vibrations 
with frequencies around 1340cm~^. Since fast interac- 
tions arise prior to the formation of the Bloch-wavc, they 
have the effect of dressing the charge with a polariza- 
tion cloud or a lattice deformation cloud. Slow interac- 
tions, on the other hand, have characteristic times much 
longer than /i/J°. They act directly on the Bloch-wave 
or the localized state. Such is the case for interactions of 
the charge carrier with low-energy intermolccular ther- 
mal phonons and librations, which in many cases can be 
considered as static with respect to the two-dimensional 
band motion ^ 50cm~^. These interactions scatter the 
Bloch-wave or localize the electronic states when the dis- 
order they introduce is large enough. An interesting dis- 
cussion of time scales can also be found in the first chap- 
ter of the book by Silinsh and Capck 

Because they dress the charge with a polarization cloud 
or lattice deformation, fast processes lead to a renormal- 
ization of the bare transfer integrals J'' and J\ and con- 
sequently increase the effective mass along all crystal di- 
rections. The case involving electron-phonon interactions 
has been discussed by several authors including Appel^ 
and Davydov4^ The purely electronic effects were treated 
in earlier work a^'^"^ in which wc calculated the renormal- 
ization effect due to the electronic polarizability in the 
bulk of the organic semiconductor. 

APPENDIX B: SPATIAL CORRELATIONS 

Consider 2 sites n and m. The polarization en- 
ergy Ep{n) depends on the spatial configuration of the 
molecules that corresponds to the polarization cloud 
around the site n. Thus, if the polarization clouds around 
the given sites n and m overlap, Ep(m) and Ep{n) are 
correlated. 

In order to evaluate the magnitude of these correla- 
tions, we calculated 

<Jn^m = {Ep{m)Ep{n)) - {Ep{m)){Ep{n)) (Bl) 

for m = n + 1 (next neighbour of n) to m = n + 11 in 
the direction of highest transfer integral. The discrete 
cluster radius was set to 18 A and the angular disorder 
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to 5 degrees. The results of these discrete ealculations are 
depicted in Fig. [5l The value for m = corresponds to 
the variance of the distribution of the polarization energy. 
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FIG. 5: Correlation between -E'p(O) and Ep{m) versus the 
number of sites between site and site m. The correlation 
decreases very fast so that we can neglect spatial correlations. 

One can readily see from these results that the correla- 



tions decrease very fast with distance to a value smaller 
than 10% of the variance. This decrease is to be re- 
lated to the induced character of the dipoles : permanent 
dipoles would lead to long-range correlations character- 
istic of Coulomb interactions. Therefore, as mentioned 
previously, we can reasonably neglect the spatial corre- 
lations of the polarization energy distribution and only 
take into account energetic correlations in the study of 
the Anderson hamiltonian of Eq. [5l 



APPENDIX C: MATRICES M„ AND M„ 

The transfer matrix equation (Eq. [9]) depends on the 
2AI X 2M matrices M„ and M„ defined as, 



n-1,0 

jn+1,0 
^n,0 



j: 



n — l.m 



n+l,m 



J, 



jn+l,M 



I E - ,5„,o 

Tn+1,0 



-1,0 



'n,0 



T-n-Hl.O 



E-S„ 



-1,1 



M„ = 



jn.m- 



7-71+ l,m 
'Jtl.ttl 



Tn+l,m 



jn,lVl —A 

Tn+l,Af-l 
^n,M-l 



n+l,M-l 



n,A/-l 



j: 



n,M-l 
n+l,M-l 



n,M-l 



V 



Tn,M+l 
'n,M 



■n+\,M 



n,M 



E — 6n,M 
jn+l.M 



rn,M+l 

'n,M 

jn+l.M 



11 



ACKNOWLEDGMENTS heavily for our calculations. They also acknowledge 

S. J. Konezny and D. B. Romero for scientific discus- 
The authors want to thank the maintainors of the sions and the Swiss National Science Fundation for the 
cluster Pleiades (pleiades . epf 1 . ch) that was used financial support (project number 200020-113254). 



* Author to whom the correspondance should be adressed, 
e-mail :libero.zuppiroli@epfl.ch 

^ P. Gosar and S.-I. Choi, Physical Review 150, 529 (1966). 
^ M.-N. Bussac, J.-D. Picon, and L. Zuppiroli, Europhysics 

Letters 66, 3 (2004). 
^ A. Troisi and G. Orlandi, J. Phys. Chem. A 110, 4065 

(2006). 

* A. A. Kostin, V. V. Savelev, and A. V. Vannikov, Phys. 
Stat. Sol. B 87, 255 (1978). 

^ N. Karl, K. H. Kraft, J. Marktanner, M. Munch, F. Schatz, 
R. Stehle, and H.-M. Unde, J. Vac. Si. Technol. A 17, 2318 
(1999). 

® W. Warta and N. Karl, Phys. Rev. B 32, 1172 (1985). 

N. Karl, Synthetic Metals 133-134, 649 (2003). 
® V. Podzorov, S. E. Sysoev, E. Loginova, V. M. Pudalov, 

and M. E. Gershenson, Appl. Phys. Lett. 83, 3504 (2003). 
^ A. Stassen, R. W. I. de Boer, N. N. losad, and A. F. Mor- 

purgo, Appl. Phys. Lett 85, 3899 (2004). 
^° N. E. Gruhn, D. A. da Silva Filho, T. G. BiU, M. Malagoh, 

V. Coropceanu, A. Kahn, and J.-L. Bredas, J. Am. Chem. 

Soc. 124, 7918 (2002). 
" V. Coropceanu, M. Malagoh, D. A. da Silva Filho, N. E. 

Gruhn, T. G. BiU, and J.-L. Bredas, Phys. Rev. Lett. 89, 

275503 (2002). 

D. A. da Silva Filho, E.-G. Kim, and J.-L. Bredas, Ad- 
vanced Materials 17, 1072 (2005). 

M. Mas-Torrent, P. Hadley, S. T. Bromley X. Ribas, 
J. Torres, M. Mas, E. Molins, J. Veciana, and C. Rovha, 
J. Am. Chem. Soc. 126, 8546 (2004). 

V. Kenkre, J. Andersen, D. Dunlap, and C. Duke, Phys. 
Rev Lett. 62, 1165 (1989). 

K. Hannewald and P. Bobbert, Phys. Rev. B 69, 075212 
(2004). 

^® K. Hannewald, V. Stojanovic, J. Schellekens, and P. Bob- 
bert, Phys. Rev. B 69, 075211 (2004). 

" A. Troisi and G. Orlandi, Phys. Rev. Lett. 96, 086601 
(2006). 

A. A. Johansson and S. Stafstrom, Phys. Rev. B 69, 235205 
(2004). 

^® L Natkaniec, E. L. Bokhenkov, B. Dorner, J. Kalus, G. A. 

Mackenzie, G. S. Pawley, V. Schnelzer, and E. F. Sheka, 

J. Phys. C : Sol. St. Phys 13, 4265 (1980). 
^° D. W. J. Cruickshank, Review of Modern Physics 30, 163 

(1958). 

21 Y. Cheng, R. Silbey D. da Silva Filho, J. Calbert, 
J. Cornil, and J. Bredas, J. Chem. Phys. 118, 8 (2003). 

22 H. M. Pastawski, L. E. F. Foa Torres, and E. Medina, 
Chem. Phys. 281, 257 (2002). 

2^ H. Houili, J.-D. Picon, M.-N. Bussac, and L. Zuppiroli, 



Appl. Phys. 100, 023702 (2006). 
2* M. Pope and C. E. Swenberg, Electronic Processes m Or- 
ganic Crystals and Polymers (Oxford Science Publication, 
1999), p. 553. 

2=5 J. K. Song, S. Y. Han, I. Chu, J. Kim, S. K. Kim, S. A. Lya- 
pustina, S. Xu, J. M. Stilles, and K. H. Bowen, J. Chem. 
Phys. 116, 4477 (2002). 

2^ M. Unge and S. Stafstrom, Synthetic Metals 139, 239 
(2003). 

2^" B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56, 1469 
(1993). 

2^ M. Unge and S. Stafstrom, Phys. Rev. B 74, 235403 (2006). 
2^ V. Oseledec, Trans. Moscow Math. Soc. 19, 197 (1968). 

G. Benettin and L. Galgani, Intrinsic Stochasticity in Plas- 
mas (Editions de Physique, Orsay, 1979), p. 93. 

^1 A. MacKinnon and B. Kramer, Phys. Rev. Lett. 47, 21 
(1981). 

^2 D. Emin, Properties of Amorphous Semiconductors (P.G. 
Le Comber and J. Mort, Academic Press, 1973), pp. 291- 
293. 

H. M. Pastawski, L. E. F. Foa Torres, and E. Medina, in^, 
p. 257. 

L N. Hulea, S. Fratini, H. Xie, C. F. Mulder, N. N. Tossad, 

G. Rastelli, S. Ciuchi, and A. F. Morpurgo, Nature Mate- 
rials 5, 982 (2006). 

V. Podzorov, E. Menard, J. A. Rogers, and M. E. Ger- 
shenson, Phys. Rev. Lett. 95, 226601 (2005). 
N. Karl, J. Marktanner, R. Stehle, and W. Warta, Syn- 
thetic Metals (1991). 

M. Kotani, K. Kakinuma, , M. Yoshimura, K. Ishii, S. Ya- 
mazaki, T. Kobori, H. Hokuyama, H. Kobayshi, and 

H. Tada, Chemical Physics 325, 160 (2006). 

R. W. L de Boer, M. Jochemsen, T. M. Klapwijk, and 
A. F. Morpurgo, J. Appl. Phys. 95, 1196 (2004). 
^® K. M. Rosso and M. Dupuis, Journal of Chemical Physics 
120, 7050 (2004). 

E. F. Valeev, V. Coropceanu, D. A. da Silva Filho, 
S. Salman, and J. L. Bredas, J. Am. Chem. Soc. 128, 9882 
(2006). 

'^^ A. von Miihlenen, N. Errien, M. Schaer, M.-N. Bussac, and 

L. Zuppiroli, unpublished (2007). 
''2 Y. Toyozawa, Prog. Theor. Phys. 12, 421 (1954). 

E. A. Silinsh and V. Capek, Organic Molecular Crystals 

(AIP Press, New York, 1994). 

V. Coropceanu and J. L. Bredas, Nature Materials 5, 929 
(2006). 

J. Appel, Solid Stat. Phys. 25, 193 (1968). 

A. S. Davydov, Theorie du solide (Edhions Mir, 1980). 



